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Abstract. We study an inverse scattering problem for the discrete Schrodingcr 
operator on the square lattice Z d , d > 2, with compactly supported potential. 
We show that the potential is uniquely reconstructed from a scattering matrix 
for a fixed energy. 

1. Introduction 

1.1. Inverse scattering. Let Z d = {n = (m, • • • 7 n d ) ; ri; e Z, 1 < i < d} be the 
square lattice, and ei = (1, 0, • • • , 0), • • • , = (0, • • • , 0, 1) the standard basis of 
Z d . Throughout the paper, we shall assume that d > 2. The Schrodinger operator 
H on Z d is defined by 

H = H + V, 

where for / = {/(n)}„ eZ d G £ 2 (Z d ) and n £ Z d 

(H f)(n) = -\^{T(n + ej) + f(n e s )} + p(n), 

3 = 1 

(Vf)(n) = V(n)f(n). 
We impose the following assumption on V: 

(A) V is real-valued, and V(n) = except for a finite number of n. 

Under this assumption, cr(H ) — a ess (H) — [0,d], and the wave operators 

(1.1) W {±) =s-lime ttR e- uRo (in£ 2 {Z d )) 

t— 5-±00 

exist and are asymptotically complete, i.e. their ranges coincide with % ac {H), the 
absolutely continuous subspace for H. Hence the scattering operator 

(1.2) S=({?<+))*?(-) 

is unitary. Associated with H , we have a unitary spectral representation 
Jb : f(Z d ) -> L 2 ((0, d);L 2 (M x );d\), 

where 

d 

(1.3) M x = |.t £ T d ; d- ^2cosx 3 ; = 2A|, 

j=i 

(1.4) T d = R d /(2nZ) d = [-tt, Tr] d . 
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Then FqS^o)* has the following direct integral representation 
(1.5) FoS(Fo)* = [ ®S{X)dX. 



o 



Here 5(A) is a unitary operator on L 2 (M\), and is called the S-matrix. 

Our main concern in this paper is the inverse scattering, i.e. reconstruction of the 
potential V from the knowledge of the S-matrix. In [10] (see also [6]), it has been 
proven that given S(X) for all energy A £ (0,d) \ Z, one can uniquely reconstruct 
the potential. 

It is worthwhile to recall the case of the continuous model, i.e. the Schrddinger 
operator —A + V(x) in L 2 (R d ). In this case, it is known that only one arbitrarily 
fixed energy A > is sufficient to reconstruct the compactly supported (and also 
exponentially decaying) potential V(x) from the S-matrix 5(A). This was proved 
for d > 3 in 1980's by Sylvester-Uhlmann [2D], Nachman [TS], Khenkin-Novikov 
[12] . There are two methods. One way is applicable to the compactly supported 
potential and based on the equivalence of the S-matrix and the Dirichlet-Neumann 
map (called D-N map hereafter) for the boundary value problem in a bounded do- 
main. The other way relies on Faddeev's theory for the multi-dimensional inverse 
scattering, in particular, on Feddeev's scattering amplitude, and allows exponen- 
tially decaying potentials. In both cases, Sylvester-Uhlmann's complex geometrical 
optics solutions to the Schrddinger equation, or Faddeev's exponentially growing 
Green function played a crucial role. (See e.g. an expositiory article [S].) However, 
since both of these methods use the complex Born approximation, the case d = 2 
remained open rather long time. Note that for the potential of the form coming 
from ellectric conductivities, the 2-dim. inverse scattering problem for a fixed en- 
ergy was solved by Nachman [IB]. See also [S]. Recently Bukhgeim [2] proved 
that, based on Carleman estimates, the D-N map determines the potential for the 
2-dim. boundary value problem. For the partial data problem, see [7]- This result 
can be applied to the inverse scattering and to derive an affirmative answer to the 
uniqueness of the potential for given potential of fixed energy. 

1.2. Main result. To study the inverse scattering from a fixed energy for the 
discrete model, we adopt the above-mentioned former approach. Namely, we as- 
sume that the potential is compactly supported, and derive the equivalence of the 
S-matrix and the D-N map in a bounded domain. 
We need to restrcit the energy in some interval. Let 



The following theorem is our main aim. 

Theorem 1.1. Fix A G I4 arbitrarily. Then from the S-matrix S(X), one can 
uniquely reconstruct the potential V . 

Our proof not only states the uniqueness, but also explains the procedure of the 
reconstuction of the potential. 

1.3. The plan of the proof. After the preparation of basic spectral results in §2 
and §3, the first task is to relate the S-matrix with the far-field pattern at infinity of 
the generalized eigenfunction of H. This is done in §4 by observing the asymptotic 
expansion at infinity of the Green operator of H . In §5, we introduce the radiation 



(1.6) 




for d = 2, 
1/2, d), for d>3. 
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condition for the Hclmholtz equation and prove the uniqueness theorem for the 
solution. We then study the spectral theory for the exterior problem in §6, with 
the aid of which we obtain in §7 the equivalence of the S-matrix and the D-N 
map for a boundary value problem in a bounded domain. The potential is then 
reconstructed from the D-N map in §8 via a constructive procedure. 

Although the main stream of the proof is the same as the continuous case, we 
need to be careful about the difference in the case of the discrete model. The first 
one is the asymptotic expansion of the resolvent at infinity. This is based on the 



stationary phase method on the surface M\ defined by (1.3 1, which is not strictly 
convex in general. This is the reason we restrict the energy on Id- The second one, 
which is more serious, occurs when we compare the far-field patterns of solutions 
to Schrodinger equations in the whole space with those of the exterior domain. 



We need a Rellich type theorem (see Theorem 5.7) and a unique continuation 
property for the discrete Helmholtz equation, which do not seem to be well-known. 
However, the former's precursor has been given by Shaban-Vainberg [19] . and the 
latter follows rather easily from it. As a byproduct, it proves the non-existence of 
embedded eigenvalues for H ([H]). We then go into the final step of computing 
the potential from the D-N map. In the continuous case, this is an elliptic Cauchy 
problem from the boundary, hence is ill-posed. However, in the discrete case, this 
is a finite dimensional problem, therefore a finite computational procedure. The 
whole proof does not depend on the space dimension. In contrast, it took a long 
time to get the 2-dim. result in the continuous case. 

1.4. Remarks for references. There are important precursors of this paper. The 
work of Eskina [6] have already announced the result of the inverse scattering for 
discrete Schrodinger operators. In particular, this paper stresses the effectiveness of 
several complex variables in the study of discrete Schrodinger operators. Shaban- 
Vainberg [19] studied the spectral theory of discrete Schrodinger operators. They 
introduced the radiation condition, proved the limiting absorption principle, and 
derived the asymptotic expansion of the resolvent at infinity including the case of 
non-convex surface. 

The computation of the D-N map for the discrete interior boundary value prob- 
lem was done in the work of Oberlin p~7|. See also Curtis-Morrow [3] and Curtis- 
Mooers-Morrow [4]. 

1.5. Notation. C's denote various constants. For any R d , x ■ y = x^yi + 
■ • • + x dVd denotes the ordinary scalar product in the Euclidean space where Xj and 
yj are j-th component of x and y respectively. For any x € R d , |a;| = (x ■ a;) 1 / 2 
is the Euclidean norm. Note that even for n — (rix, ■ ■ ■ , rid) € Z d , we use |n| = 

(Y^rj=i l^'l 2 ) 1 ^ 2 - For two Banach spaces X and Y, B(A, Y) denotes the totality 
of bounded operators from X to Y. For a self-adjoint operator A on a Hilbert 
space, (J (A), cr ess (A), Udi SC {A), <r ac (A) and <J P {A) denote its spectrum, essential 
spectrum, discrete spectrum, absolutely continuous spectrum and point spectrum, 
respectively. For a set S, denotes the number of elements in S. We use the 
notation 

(t) = (l + i 2 ) 1/2 , teR. 



1.6. Acknowledgement. The authors are indebted to Evgeny Korotyaev for use- 
ful discussions and encouragements. The second author is supported by the Japan 
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Society for the Promotion of Science under the Grant-in- Aid for Research Fellow 
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2. Momentum representation 

2.1. Discrete Fouri er t ransform. From the view point of dynamics on the lat- 
tice, the torus T d in ( 1.4) plays the role of momentum space. Let U be the unitary 
operator from £ 2 (Z d ) to L 2 (T d ) defined by 

(W f)(x) = (27r)- d / 2 J2 f(n)e- in -*. 

n£Z d 

Using this discrete Fourier transformation, the Hamiltonian H is represented by 

H = UHU*=H a + V, H =UH U* 1 V =UVU*, 
where Hq is the multiplcation operator: 

d 



(2.1) H = -(d-^cosxj) =: h(x), 

and V is the convolution operator 

(Vu)(x) = (2tt)-^ 2 f V(x - y)u{y)dy, V(x) = (2tt)-^ 2 V V{n) e - in - x . 

2.2. Sobolev and Besov spaces. We define operators Nj and Nj by 

(Nj /) (n) = nj f(n) , Nj = UN U* = i A . 
We put iV = (iVi, • • ■ , Nd), and let iV 2 be the self-adjont operator defined by 

d 

N 2 = N 2 = -A, on T d , 

j=i 

where A denotes the Laplacian on T d = [— ir, ir] d with periodic boundary condition. 
We put 

\N\ = VN 2 = \/ z A. 

For s € R, let W be the completion of D(\N\ S ) with respect to the norm \\u\\ s = 
\\(Nyu\\ : 

H s = {ue V'(T d ) ; = ||(A0 S U || < oo}, 
where V'(T d ) denotes the space of distribution on T d . Put U = U a = L 2 (T d ). 

For a self-adjoint operator T, let x( a < T < b) denote the operator Xi(T), 
where Xi(ty is the characteristic function of the interval / = [a, b). The operators 
x(T < a) and x(T > b) are defined similarly. Using the series {rj}° c L with r_i = 0, 

= 2 J (j > 0), we define the Besov space B by 

oo 

B = {/ G H; \\f\\ B = E r i /2 Hx(^-i < |JV| < r,-)/|| < oo}. 

Its dual space B* is the completion of % by the following norm 
IMIb* = s\iprT 1/2 \\ x (r j - 1 < \N\ < rj )u\\. 

j>0 
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The following Lemma 2.1 is proved in the same way as in [TJ. 

Lemma 2.1. (1) There exists a constant C > such that 

( 1 \ 1/2 

C- x \\u\\ B . < f sup -||x(|JV| < R)ufj < C\\u l3 . 

Therefore, in the following, we use 

( 1 -W2 

Mb* = sup-||x(|iV| <i?)«| 



\R>\ R 
as a norm on B* . 

(2) For s > 1/2, the following inclusion relations hold : 

H s c B c H 1/2 cHc -H~ 1/2 cB* cH 



We also put % = ^ 2 (Z d ), and define H s , B, B* by replacing N by AT. Note that 
T~L S = U*H S and so on. In particular, Parseval's formula implies that 



u\\ ns = \\u 



2 -^%= £(i + M 2 ) s R-)l 2 



nez d 

2 



» = sup 4 X! i^' 1 



1^*11 K ■ ;-) 

|n|<_R 



w(n) being the Fourier coefficient of u(x). 

2.3. Resolvent estimate. 

Lemma 2.2. (1) a(H Q ) = o- ac (H ) = [0,d]. 
cr eS5 ( J ff) = [0,d] ! a fa (fl)cR\[0,d], 

(s; o- p (fl)n((o,d)\z) = 0. 



Proof. The assertions (1), (2) follow from (2.1 1 and Weyl's theorem. The asser- 
tion (3) is proven in |llj . □ 

Let R(z) = (H - z)' 1 . 

Theorem 2.3. (1) Let s > 1/2 and A 6 (0, d) \ Z. TTien i/iere exists a norm limit 
R{\ ± iff) := lim e ->o #(A ± ie) € B("H s ;-fr s ). Moreover, we /iawe 

(2.2) sup||i?(A±iO)|| B(g .^ <oo. 

for any compact interval J in (0, d) \ Z. TTie mapping (0, rf) \ Z 9 A 4 i?(A ± iO) 
«s norm continuous in B(H S ;"H _S ) and weakly continuous in B(£? ;£?*). 
(2) H has no singular continuous spectrum. 

For the proof of Theorem |2.3| see Lemma 2.5 and Theorem 2.6 of [10]. Note that 

(2.3) Vfc(ar) = h(x) e {0, 1, ■ ■ • ,d}. 

This is the reason why the set of thresholds {0, 1, • • • , d} appears. 

3. Spectral representations and S-matrices 
We recall spectral representations and S-matrices derived in §3 of [TO]. 
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3.1. Spectral representation on the torus. We begin with the spectral repre- 
sentation in the momentum space. Let us note 

i ^ ^ 

K x ) = g ( d ~ cos x i) = sin2 (y) > 

which suggests that the variables y = (yi, ■ ■ ■ ,yd) £ [— 1, 

2/j = sin — , Xj = 2 arcsin g/j 
are convenient to describe Ho. Note that for A € (0, d) \ Z 

(3.1) x(-\/A0) = (2 arcsin(yA6» 1 ), • • • , 2 arcsin (VA0 d )) , G 

gives a parametric representation of 

(3.2) M A = {ieT d ;/i(i)=A}. 
We equip M\ with the measure 

dM A = J(y f X9)d9, 

x(y) being the characteristic function of [— 1, l] d . Then we have 

rfa; = J{y)dy = dM^dX, dM\ = |V _^ A , 

where dM>, is the measure on M\ induced from dx. Let L 2 (M\) be the Hilbcrt 
space with inner product 



(^■0)l 2 (m a ) = / (pipdM x . 
Jm x 

We define Jb(A)/ = Ttm a /, where TrM A is the trace on M\. More precisely, 
(3.3) (-Fo(A)/) (^) = /(s(VA0)). 

It then follows for Rq(z) — (Hq — z)^ 1 

^-((Ro(\ + i0)-Ro(\-i0))f,g) L 2 iT «) = (FoWfMXMwiMx), 



for A G (0, d) \ Z and /,g G C 1 ^). We then have by Q 
(3.4) J„(A)eB(6;L 2 (M A )). 

Using this formula, we can derive the spectral representations of i?o an d -ff ■ How- 
ever, we omit it. 



DISCRETE SCHRODINGER OPERATORS 



7 



3.2. Spectral representation on the lattice. Wc define the distribution S(h(x) 
A) e V'{T d ) by 



f(x)S(h(x) - X)dx 



f(x)dM x , /GC°°(T^ 



Ms 



Then, from the definition of J-q(X)*: 

(.Fo(A)/, 0) l2(Ma) = (/, J"o(A)*0) i2(Td) , 

we see that J-o(A)* defines a distribution on T d by the following formula 

T {X)*(j) = 4>(x)8(h(x) - A). 

Here the right-hand side makes sense when, for example, <fi g C°°(M\) and is 
extended to a C°°-function near M A . Then J-o(A)*(/> = U* <f) is computed as 



(2tt) 



d/2 / Jn-x 



(3.5) 



= (27T)- d / 2 

=(2^)- d / 2 



-pel 



(x)S(h(x) — X)dx 



(x)dM x 



S d-i 2 



-J(VX9)d9. 



In the lattice space, we define V (0) (A,6») = {V> (0) (™, A, 6>)} n6Zd , where 



^ (0) (n,A,6») = (2ir)- d/: 



(Va) 



d-2 



D m-a;( y/XO) 



J{VX9) 



(3.6) 



(27r)- d / 2 2 d - 1 (\/A) d - 2 x(VA6l) 



3 m-:c(\/A0) 



Here x(y) is the characteristic function of [—1, and x(\f\Q) is defined by (3.1 ). 
By ([3j3} and we have for <j> e L 2 {M\) 

(^ O (A)*0)(») = (2rr)- d / 2 / e* 1 -^) dM A 

^ (o) (^,A,6»)0(x(v / A6»))d6i. 

-1 

We can also see for rapidly decreasing / on Z d 

nGZ d 

The spectral representation for H is constructed as follows. We put 
(3.7) = ^o(A) (l-VR(\±iO)) , Ae(0,d)\Z. 



Then by and ( |2.2| ) 

^( ± )(A)eB(S;L 2 (Af A )). 
We define the operator J^) by (.F (±) /)(A) = -F (±) (A)/ for / eg. 
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Theorem 3.1. (1) is uniquely extended to a partial isometry with initial set 

H ac (H) and final set L 2 (T d ). Moreover it diagonalizes H: 

(3.8) {^Hf){\) = \(^f){X), f€H ac (H). 

(2) The following inversion formula holds: 

(3.9) 7=s-lim/ T {± \X)*(T^f)(\)d\, feH ac (H), 

JV->oo J In 

where In is a union of compact intervals in (0, d) \ Z such that In — > (0, d) \ Z. 

(3) T (± \\)* E B(L 2 (M A ) ; B*) is an eigenoperator for H in the sense that 

{H -\)F {±) {\)*<P = Q, <t>^L 2 {M x ). 

(4) The wave operators 

W {±) =s-\ime ItS e- ItSo 

t->±oo 

exist and are complete. Moreover, 

= (-F (±) )*^o. 

3.3. Scattering matrix. The scattering operator S is defined by 

S= (W + )*W— 
We conjugate it by the spectral representation. Let 

S = FoS^o)* , 

which is unitary on L 2 ((0, d); L 2 (M\);dX). Since S commutes with H , S is written 
as a direct integral 



/ 

J(o 



}S(X)d\. 

'(0,d) 

The S-matrix, 5(A), is unitary on L 2 {M\) and has the following representation. 
Theorem 3.2. Let A e (0, d) \ Z. Then 5(A) is written as 

S{\) = 1 - 2ttL4(A), 

where 

(3.10) A(X) = Jo(A) (l-VR(\ + iQ)) KF (A)* = ^ (+) (A)^o(A)*, 

and is called the scattering amplitude. 

4. Asymptotic expansion of the resolvent at infinity 

4.1. Stationary phase method on a surface. Let S be a compact C°°-surfacc 
in R d of codimension 1, and dS the measure on S induced from the Euclidean 
metric. For a(x) € C°°(S) and k e R d , we put 

(4.1) J(jfe) = / e M - fc a(x)dS*. 
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Theorem 4.1. Let N(x) be an outward unit normal field on S, and W{x), K(x) 
the Weingarten map and the Gaussian curvature at x G S, respectively. Assume 
that there exists a finite number of points x± G S, j = 1, • • ■ , v , such that 



k/\k\=±N(x ( £ ) ), 



ano 



d that K(x^£ > ) =/= 0. j = 1, • • ■ , v. Then we have as p = \k\ — > oo 



(4.2) 



I{k) = p-^' 2 ^^"^ A+{xf) 

3=1 

V 

, n -(d-l)/2 ST p *k-x ( J } A ( Jj)s , n( ( d +l)/2s 



-D/2 £ ert-^^a.W)) + 0(p-(rf+D/2) ) 

3=1 

w/iere 

(4.3) At (a) - (27r)(' i - 1 )/ 2 |X(2;)|- 1 / 2 e =FssnW ^ ,ri / 4 a(a;). 

and sgnVF(x) = n + — n_, n + feeing i/ie number of positive (negative) eigen- 

values of W(x). 

For the proof, see Lemma 2.2 and appendix of Q3]. See also [13]. If S is 
represented by Xd = f(x'), x' = (xi, ■ ■ ■ ,Xd-i), the Gaussian curvature is given by 

'' ' df f _,^ 3 , ^-(<i+i)/2 ^ ^ 
For d = 2, the Gaussian curvature of the curve f(xi, x%) = is computed as 

\fx2X2 ' /xi — 2fx 1 x 2 ' fxifxz + fxixi ' /i 2 | 



,4.4, KW _(£(V M) V.)-^--(- * M 



(4.5) ^(x 1; x 2 ) 



(fl+fl) 3/2 



4.2. Convexity of M\. As will be seen below, the shape of depends highly 
on the space dimension and A. We know that V/i(x) ^ on M\ if A g" Z. Assume 
that at a point in M>, dh/dxd = (sinxd)/2 ^ 0. We take Xi, ■ ■ ■ ,Xd-i as local 
coordinates, and differentiate h[x) — A to get 

dxd n 

<9x,j 3x rf . <9 2 x d 

djj cosxj + cosx d — - h sinx d - — - — = 0, 

oxi axj oxiOXj 

for i, j — 1, ■ ■ ■ , d — 1. We put <p = ^ ■_ 1 fcjX.,-. Then we have on M\ 

dip dxd smXi 

kj -\- kd. — kj kd. 



oxi oxi sin 

, 9 2 a; d k d , , , , 2 , • • \ 

fcd — = — -r^ rg (Oij COS Xj (SUl Xd) + Sin Xi Sin Xj COS Xd j • 



dxidxj dxidxj (sinx^) 3 
Suppose dip/dxi = 0, i = 1, ■ ■ ■ ,d—l. Then 

ki = p sin Xi, i = 1 , • • • , d, 
p = \k\ ((sinxO 2 + • • • + (sinx d ) 2 ) _1/2 
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Therefore we have 



(4.6) 



d 2 ip _ 1 
dxidxj (sinxd) 2 /? 



(Sijkj cosxj + kikj cosxd) 



Now let us compute the determinant det (d 2 <p/ dxidxj) . 

(1) The case d — 2. Using ki — psinxi, we have 

k\ cosxi + k\ cos 2:2 = (?{cosx\ + cosa^Xl — cosxi cos 2:2) 

= 2p 2 (l — A)(l — cosxi cosa^). 

Since A ^ 1, this vanishes if and only if cosxi = cos X2 = ±1, i.e. x\ = or n, and 
X2 = or 7T. However in this case, h(x) = sin 2 (cci/2) € Z. This implies that 

(4.7) d 2 (p/dxj ^ for A e (0, 1) U (1, 2). 

Therefore M\ is a closed curve in T 2 , and convex in the fundamental domain 
R 2 /(27rZ) 2 , as is seen from the figures (Figures 1, 2, 3) below. Let us remark here, 
in view of Figure 3, in the case 1 < A < 2, it is convenient to shift the fundamental 
domain so that R 2 /(2ttZ) 2 = [0,2vr] 2 . 





-3 



-10 1 2 3 



-3 



-1 1 



Figure 1. d 
2, A = 0.25. 



Figure 2. d 
2, A = 0.75. 




-6 -+-2 2 4 6 



Figure 3. d = 2, A = 1.25. 
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(2) The case d = 3. By a direct computation, wc have 
det (dijk\ cosxj + fe,-fcj cosa^ 

= (k\ cos Xi cos x-i + /cf cos X3 cos Xi + cos X\ cos X2) , 

which can vanish when e.g. cosa^i = cosa;2 = 0, cos x$ — 1/2. Therefore in 3- 
dimensions, M\ may not be convex. The following Figures 4, 5, 6 explain the 
situation in 3-dimcnsions. 




Figure 4. d 
3,A = 0.45. 




Figure 5. d = 
3,A = 2.55. 





Figure 6. d = 3, A = 1.45. 



Here, we note the following simple lemma. 

Lemma 4.2. If —1 < j/j < 1, i = 1, • • • , d, and d — 1 < y\ + • • • + yd < d, we have 
y t >0,i = l,---d. 



Proof. Suppose e.g. yd < 0. Then 

Vi + V2 H h Vd < Vi + 

which is a contradiction. 



Vd-i <d-l, 



a 
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By (|4.6[), we have 



d-i 



4-^ dxidx-j 

,.7 = 1 J 



'J 



^ / d-1 d-1 2 \ 



which has a definite sign if cosxi > 0, i = 1, • • • , d and sin aid > 0. By virtue of 
Lemma 4.2, it happens for < A < 1/2. Let us also note that for d — 1/2 < A < d, 
we have the same conclusion since cosxi < (i = 1, • • • ,d), sinxd < 0. Recall 
that when d > 3 the definition of the Gaussian curvature depends on the choice 
of direction of the unit normal N(x) on S. We choose N(x) in such a way that 
K(x) > on S. 

With this convention, we have proven the following lemma. Recall the interval 



Id defined by ( 1.6 ) 



Lemma 4.3. If A £ Id, all the principal curvature of M\ are positive. 

As has been noted above, in the case 1 < A < 2 (d = 2) or d— 1/2 < A < d (d > 3), 
we should shift the fundamental domain so that R d /(27rZ) d = [0, 27r] d (See Figures 
3, 4, 5). To fix the idea, in the sequel, we deal with the case R d /(2nZ) d = T d = 
[-n,ir] d . 

Under the assumption of Lemma 4.3, M\ is strictly convex. Let N(x) be the 
unit normal field on M\ specified as above. Then for any to £ S^ 1 , there exists a 
unique pair of points x±(X,lu) in M\ such that 

(4.8) N(x±(\,u))) = ±lj. 

Since N(—x) — —N(x), we see that X-(\,w) — —x+(\,cj). Therefore, we let 

(4.9) x±(\,w) = ±x 00 (X,lj). 

We can now compute the asymptotic expansion of the free resolvent 

(4.10) (Ro(z)f) (m) - r °( m - n ' 

n£Z d 



(4.11) r (k,z) = (2TT)- d — dx. 

J T d n(x) — z 

We put 

(4.12) w k = k/\k\, keR d \{0}. 
Lemma 4.4. Assume X G Ij. Then we have as \k\ — > oo 

r o (k,X±i0) 

= ± l ( 27r |fc|)-(^i)/2 e ± l (' £ -^(^^)-(' i - 1 )-/4) KixiiX,^))- 1 ^ + om -(d+D/2y 

\v x h(x±(X,u k ))\ 



Proof. Take e > small enough so that 

'(0,1), d=2, 



(A — 2e, A + 2e) C 



(0,1/2), d>3. 



Let x(t) G Co°(R) be such that x(t) = 1 for \t\ < e/2, x(t) = for \t\ > e, and 
assume that |Re z — A| < e/4. We split r (fc, z) into two parts 

r (k,z) = A(k,z) + B(k,z), 
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A(k,z) = (2n)- d ( X ^~ X) e ik - x dx. 
J T d h(x) — z 

Then, by integration by parts, for all N > 

B(k,z)=0(\k\- N ), |fc|->oo. 
Letting S(t) — {x e T d ; h(x) = t}, wc write A(k,z) as 

A(k,z) = (27r)" d f X+e fl^dt, a(t,k) = I e ^^-\dS(t). 



(It. 



We then have 

(4.13) _ j _^ < ft = ± Mr0 (A > *)+p.v.y^ — 

By Theorem 4.1, for f G (A — e, A + e), a(t, k) admits the asymptotic expansion 
a(t,k) -a (t,fc) + O(|fcr (d+1)/2 ), 



V|fc|/ |V 2; /i(a; + (t,a; fe ))| 
(4-14) ^/27r\W-i)/a ifc . aM( )+((J _ 1)wi/4 ^K(x-{t, Uk ))-V 2 



|*|/ |V x MaJ-(*.Wfc))| 

where x±(i,Wfc) is a stationary phase point on S(t). 

We compute the asymptotic expansion of the 2nd term of the right-hand side of 
(4.13). Differentiating h{x±(t, u> k )) = *> we have 

V x h(x±(t,Uk)) ■ d t x±(t,oj k ) = 1. 

Therefore, letting 

S =UJ k ■ X ± (t,U> k ) - LO k ■ X±(A,Wfe), 

we have 

ds V x h(x ± (t,ui k )) 1 

— = uj k ■ d t x ± {t,uj k ) = — — T ■ d t x±{t,u k ) = — — T ^j, 

at \\/ x h(x±(t,uj k )\ \\/ x h(x±(t,L} k ))\ 

which implies 

t - A = s\V x h(x±(\,u> k ))\ + 0{s 2 ). 

We then have 

1 X (t- ^Kjx^t,^))- 1 / 2 dt = b ± (s,u k ) 
t — \ \\7 x h(x±(t,uj k )\ ds s 

where b±(s,uj k ) is a smooth function such that 

0±(U,Wfe) - — rry. 

\V x h(X±{X,UJ k ))\ 

Taking S > small enough, we have by integration by parts 

f s e ±i\k\s f \k\s ■ 
p.v. / b ± {s,u k )ds = ±2i dsb±(0,u) k ) + O(\k\- x ) 

J-5 S J S 

= ±TTib±(0,UJ k )+O(\k\- 1 ), 
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which implies 
(4.15) 



p.v. 



A-e t — X 

W\ 



dt 



e ±ife-a; 00 (A,aj fc )=F(<i-l)^/4p y 



S e ±i\k\s 
■S s 



b ± (s,u k )ds + 0(\k\^ d+1 ^ 2 ) 



\\k\J *~ '■' ' N 



\Vxh(x±(\, u) k ))\ 



Plugging (4.131, (4.141 and (4.151, we obtain the lemma 



□ 



Lemma 4.5. We have as \m\ —> oo 

(m- n) ■ x±(X,(jj m - n ) = (m -n) ■ x±(\,u} m ) + 0(|m| _1 ). 

Proof. We extend x±(A, fc) as a function of homogeneous degree in k. Letting 
e = l/\m\, we have 

u m - n = (w m - era)/|u; m - en| = u m + e((w m ■ n)w m - n) + 0(e 2 )- 
Using /i(j;±(A, u; m _„)) = A, we have 

V x /i(a;±(A,a; m _ n )) • ^-x± (A, w m _„) = 0. 
Since V x /i(x±(A, w)) is parallel to w, we then have 

u m ■ —X±(X,CJ m -n) = 0, 
ae £=o 

which implies 

m ■ a;±(A,w m _„) = m ■ x±(X,u m ) + 0(|m| _1 ), 
and the lemma follows immediately. □ 

Lemmas 4.4 and 4.5 imply the following lemma. 
Lemma 4.6. If X E Id and f{n) is compactly supported, we have as \k\ — > oo 
(i?o(A±zO)j) (k) 

n 

+ o(|fc|-( d+1 )/ 2 ), 

K(x±(X,tjj k ))~ 1 / 2 



(4.16) 



a±{X,uJk) 



\V x h(x±(X,u> k ))\ 



Recalling the definition of x{\/X6) in (3.1) and the fact that the Gauss map 
is a diffcomorphism for a strictly convex surface, define 6{X,oS) by the relation 
x{\/X0(X,lo)) = Xoo(X,lj), i.e. 



(4.17) 6j (A, ui) = i sin x^j (A, j , j 



= !,••• ,d. 
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We define the reparametrized Fourier transforms C?o(A) and t/' ± -'(A) by 

(4.18) (Go(X)f) (w) = (?oMf) («(A,w)), 

(4.19) £ (±) (A) = 0o(A)(l - VR(X ± tO)). 



Lemma 4.6, the definition (3.7) and the resolvent equation imply the following 
theorem. 

Theorem 4.7. If \ E Id and f(n) is compactly supported, we have as \k\ —¥ oo 
(i?(A±*0)/) (A) 

= e ±(3 - d) " /4 ^|fc|- (d - 1)/2 e ± * a; ~ (A ^ ) a ± (A,^) (^ (±) (A)/) (±w fc ) 
+ 0(|jfe|-< d+1 )/ 2 ). 

5. Radiation conditions on Z d 



The aim of this section is to introduce the radiation condition (Definition 5.5 ) 



and prove the uniqueness theorem (Theorem 5.9 ) 



5.1. Green's formula. For m,n € Z d , we write m ~ n, if |m — n| = 1, i.e. there 
exists j such that m = n ± ej. We define the discrete Laplacian A<2i sc on Z d by 

(5.1) (Adi sc tt)(n) = -(H u)(n) = - ^ (u{m) -u(n)). 

A set D C Z rf is said to be connected if for any m, n € -D, there exist m^) € £>, 
j = 0, ■ • • , k such that m* 3 ') ~ m^ +1 ), j = 0, • • • , — 1, and = m, mS k > = n. 
A connected subset D C Z d is called a domain. For a domain Z? C Z d , we define 

(5.2) deg (n) = *{m € I? ; m ~ n} , n £ D, 

(5.3) I)= {ne L>; deg(n) =2d}, 

(5.4) &D = {n € L>; deg (n) < 2d}. 

The normal derivative at the boundary is defined by 

(5.5) (d? u) (n) = - 51 ( S ( n ) " "( m ))' n e 

o 



Note that, compared with (5.1|, m and n are interchanged. Then the following 



Green's formula holds (see e.g [5] and [XT]): 

^ ({A disc u)(n) ■ v(n) - u(ra) • (Ad isc w)(n)) 

(5.6) ne ° D 

= ((^fi)(n)-«(n)-u(n).(^«)(n)). 
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5.2. Radiation condition. For m, n such that m ~ n, wc define the difference 
operator 9 m _„ by 

(<9 m _ n /)(n) = f(m) - f(n). 
Lemma 5.1. (1) Let n(s) = n + s(m — n), where m ~ n. T/ien we /icwe 

d m -n(n ■ x x (\,L) n )) = / (m- n) • ^oo(A,w„( s ))rfs. 
Jo 

Ifm^n, we have as \n\ — > oo 

d m - n (n ■ Xoo(A,a; n )) = (m - n) • Xoo(A, w n ) + 0(|n| _1 ), 
9m-n (e in - x °° ( - X ' u «A = ^ e i(m-«)^oo(A,a;„) _ ^ e in- Xoo (A,cc„) + 0(| n |-l). 



Proof. Differentiating /ifioofA, w n j- s -j)) = A, we have 

(V x /l)0Coo(A,W n ( s ))) • ^-Xoc(X,U n ( s )) = 0. 

Since (V x /i)(a; 0O (A, w n ( s ))) is parallel to n(s), we then have 

' ^oo(A,w„ (s) ) = 0, 

which implies 

^(n(s) • Xoo(A, w„( s ))) = (m- n) ■ x 00 (X,uj n ^). 

Integrating this equality, we obtain (1). Since w„( s ) = u n + 0(|n| _1 ), (2) follows 
from (1). □ 

We now introduce the rectangular domain D(R) such that 

(5.7) D{R)= {n G Z d ; n G [-R,R} d }, R>0, 
and the radial derivative d ra d by 

(5.8) (d rad u){k) = - A (u(m)-u(fc)), 

meaZ3(i?,(fc)),m~fc 

(5.9) i?(fc) = max 1/U, k G Z d . 

1<J<<2 J 

We put 

(5.10) v4 ± (A, Wfc ) = ^ ^ ^ e ±i(m-fc).».(A^) _ ^ c fe = A. 



Lemma 5.2. fij TTie right-hand side of (5.10) does not depend on \k\. 



(2) There exists a constant e (A) > such that 

±lmA ± (X,cj k ) > e (A), 

/or any ui k . 
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Proof. If m ~ k, m E dD(R(k)), then m — k = ±ej for some j. This ±ej 
depends only on uj k} which proves (1). 

Recall that Vh(x) = g(sin:ri, • • • , sina^), hence letting ojkj be the j-th compo- 
nent of cjfc, we have 

sin(a; 00: ,-(A,cj fc )) = aj k ,j 

for some constant c > 0. Suppose m <~ k, m E dD(R(k)). If Wfej > 0, then either 
mj = kj or mj = kj + 1. If w^j < 0, then either mj = kj or mj = kj — 1. We then 
have that sin((m — fc) • x oc (X,uJk)) = c\^>k,j\ for some j such that w^.j 7^ 0. Since 

±ImA±(A,u;fc) = * ^ sin((m - fc) • Xoo(A, w fc )) 

mGaD(_R(fc)),m~fe 



I E i wfe .ji' 

medD(R(k)),m~k 



and ^ . 0J k j = l, the lemma follows. □ 
Let us introduce two auxiliary norms, £?j^-norm and fJ^-norm, on ^* by 
llwllk - sup i ^ Rn)| 2 , 



1 8* ~" " J 
° R fl>l,ij£R 

n£_D(_R) 



° z P >i, P ez P 



E l"W| 2 - 

n£D(p) 

Lemma 5.3. These three norms || • \\ ■ \\g, , and \\ ■ are equivalent. 

Proof. Let A{R) = {z E C d ; (£)J =1 \z 3 \ 2 ) 1/2 < R}, B(R) = {z E C d ; ma Xj \ Zj \ < 
R}. Then there is a constant 5 > such that A(6R) C B(R) C A(R/6), Vi? > 0. 
This implies 

± E ^ ^ E i-c-)i 2 < ^ E i^)i 2 - 

|n|<«fl n£D( ° fl) |n|<fl/« 

Taking the supremum with respect to R > S or R > 1/5, we get the equivalence of 
|| • norm and || • ||g-» norm. 

Next we show the equivalence of || • ||g, norm and || • norm. Note that 

f( r ) = ° l"( n ) | 2 is a right-continuous non-decreasing step function on (0, 00) 

n£D(r) 

with jump at integers. For R > 1, we take p(R) = [R] = the largest positive integer 
such that p(R) < R- Then we have 

SUP 7? E l^(")| 2 < SU P -7m E l"( n )| 2 ' 

a>i R r>i P\R) 

n£D(R) n£D(p(R)) 

The converse inequality is proven by the following inequality 

Sup "Tm E l"(«)| 2 <sup- V |u(n)| 2 . □ 
R>1 P{R) R>1 R o 

neD(p(R)) n£D(R) 
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Lemma 5.4. (1) If f e £°°(Z d ) satisfies |/(n)| < C(l + \n\)-( d -^ 2 , then 

(5.11) sup i V |/(n)| 2 < oo, i.e. /eB*. 

i it '—^ 



R> 



\n\<R 



(2) If |/(n)| < C(l + Inl)-^- 1 )/ 2 - 6 , e > 0, ffcera 

(5.12) lim i ^ |/H| 2 = 0. 

|n|<fl 

Proof. We compute the norm ■ We first show 

(5.13) J2 I/M| 2 = 0(1), 

neD°(p)\D(p-l) 

o o 

as p — > oo. In fact, for any p € Z, p > 1 and n ED(p) \ D(p — 1), we have 
p-1 < |n| < Since *{ n £ D°(p) \ D{p°- 1) } = {2p+ l) d - (2p- l) d < Cp d -\ 

E I/HI' ^ Gp-^- 1 ) *{n G \ 1) } < C. 

n£Z»(p)\D(p°-l) 

On the other hand, since 



E i/>)i 2 = E E 
p=i 



l/HI 2 + l/(o)l ; 



ieD(R) 



n£D(p)\D(p-l) 



for every positive integer R, we have ^ o \f(n)\ 2 = O(R) by (|5.13|). This 



proves (1) by Lemma 5.3 



Assume |/(n)| < C(l + |n|) ^ 1 - ) / 2 e for some e > 0. By the similar computa- 



tion, we have ° l/( n )| 2 = which proves (2). 

nED(R) 

For /, g € 23*, we write 
(5.14) 7~£<=>- lim i E \f(n)-g(n)\ 2 =0. 

R— >oo it ' — ' 

|n|<R 



□ 



As we have seen above, (5.14) is equivalent to 



lim i E l/W-?(n)| 2 =0. 
it— >oo it ' — ' 

Now let us consider the equation on Z d : 
(5.15) (H-X)u = f. 



Definition 5.5. A solution u(k) € S* of (5.15) is said to be outgoing (for +) or 
incoming (for — ) if it satisfies 

(5.16) (d rad u)(k) ^ A±(\,u k )ti(k), 



in the sense of (5.14) 
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Theorem 5.6. Let X € Id- If f is compactly supported, R(X±iO)f is an outgoing 
(for +) or incoming (for —) solution of the eguation (H — X)u = f. 

Proof. Since loofA,^) is homogeneous of degree in k (see also the proof of 
Lemma 4.5), we have as |fc| — > oo 

a;oo(A, Lo k±ej ) = loofA.Wfc) + 0(|fc| _1 ). 
Then we have for any fixed n G Z d 

(5.17) e ±in-x x (\, Uk±ej ) _ e ±in-x x (X, Uk ) = OQk]- 1 ), \k\->QO. 

If / is compactly supported, (G^\X)f) (u>k) is smooth with respect to k, so that 
we have from (5.17) 

(5.18) d m - k ( a± (A, Wfe )(g( ± )(A)./)(± Wfc )) = Odfcr 1 ). 

We put = R(X ± iO)f. Theorem [47] yields 

(d r a d u {±) )(k) 

(5.19) =C ± \k\-^^ ( Yl (dm-k^m) +0(\k\-^/% 
as \k\ — > 00, where 

^\k) = e ±lk - x ^ x ^a ± (X,u Jk )(^ ± \X)f)(±Lo k ). 
Lemma 5.1 (2) and ( |5.18[ ) imply the theorem. □ 

5.3. Rellich type theorem. The following is an analogue of the Rellich type 
theorem for Schrodinger operators in H d ([18 ). 

Theorem 5.7. Let X £ (0,d) \ Z. Suppose a sequence {u(n)} defined for \n\ > 
Ro > satisfies 

(-A disc - X)u = 0, |n|>i? , 
lim i J] |£(n)| 2 = 0. 

R <\n\<R 

Then there exists R\ > Rq such that u(n) = for \n\ > 
For the proof, see [TT], Theorem 1.1. 

5.4. Uniqueness theorem. 

Theorem 5.8. Let X G Id, and suppose that f is compactly supported. Let be 
the outgoing (for +) or incoming (for —) solution of the equation (H — X)u^ = f. 
Then 

(u (±) J)-(lu {±) ) = 2 i lim V ImA ± (X,oj k )\u^(k)\ 2 . 

R— >oo £ — * 

keD(R)\D(R-l) 
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Proof. By Green's formula, we have 

]T ((A dlsc u^)(k) ■ m{kj - u&(k) ■ (A disc «(±))(fe)) 



= ((^ ( " ) M( ± ))(fc)-u(±)(fc)-w( ± )(fc)-(^ (p) u(±))(fc) 

The left-hand side converges to (u^\ /) — (/, u^) by the equation. Changing the 
order of the summation, we can see that the right-hand side is equal to 



kedD(p) 



\ J2 J2 (u (±) (fc)-w( ± )(m)-M (±) M-« (±) (fc) N 

P) ° 

J2 ((a ra / ) )H-PM-fi (±) H'P< i W) 



rneD( P )\D{p-l) 

As p — > oo, we can replace draju^ by A± (A, w^m^), and prove the theorem. □ 

Th eore m 5.9. Let A G 7^. /// is compactly supported, then the outgoing solution 
of ( 5.15) is unique and given by R(X + iO)f- The incoming solution is also unique 
and given by R(X — iO)f. 



Proof. In view of Theorem |5.6[ we have only to prove t he u niqueness. Let u be 
the outgoing solution of (H — X)u = 0. Then, by Theorem 
we have lim^^oo Y] o □ |u(fc)| 2 = 0. This implies 



5.8 



and Lemma 



5.2 



(2), 



k£D(Ft.)\D(R-l) 

lim - 

B.^oo R 



lim I E l"WI 2 =0, 



keD(R) 



i.e. u — 0. We can then use the Theorem 5.7 and the unique continuation theorem 
(see [IT], Theorem 2.1) to see that u = 0. □ 

6. Exterior problem 

6.1. Helmholtz equation in an exterior domain. Let D(R) be a rectangular 
domain in (5.7), and take a sufficiently large integer i?o > such that 

(6.1) 
We put 
(6.2) 

(6.3) 



suppV CD(R ) 



^int — D(R ), 



Therefore Q int = D{R ) = [-Rq, R„] d n Z d , and 

d 

(6.4) dQ int = dn ext = \J {n; \m\ < Rq, (i^j), 1^1 = ^ + 1}. 

i=i 

The spaces B, B* and T-L s on Qext are defined in the same way as in the whole 
space. Let H ext = — A^sc on fi exi with Dirichlet boundary condition. 
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Lemma 6.1. (1) H ext is self-adjoint, and a(H ext ) = [0, d]. 

(2) opCff« rt )n((o > d)\z) = 0. 

Proof. The assertion (1) follows from the standard perturbation theory, and (2) 
is proved in Theorem 2.4 of [11]. □ 

^ o 

For the solution of the equation (— A<& sc — X)u = / in f2 e xt> the radiation con- 
dition is defined in the same way as in §5. The following theorem is proved in the 



same way as in Theorem 5.9 



Theorem 6.2. Let A £ Id- Then the solution of the equation (—Adi SC — X)u = 

o 

in £l ex t, satisfying the Dirichlet boundary contidion and the outgoing (or incoming) 
radiation condition vanishes identically on Q e xt- 

We prove the limiting absorption principle for R ext (z) — (H ext — z)^ 1 . 

Theorem 6.3. (1) For X £ Id and f £ B, the weak *-limit exists 

lim R ext {\ ± ^e)7=: R ext (X ± i0)fe B* . 

(2) For any compact set J C Id, there exists a constant C > such that 

\\R ext (\±iO)f\\ s * <c\\f\\ s , XeJ. 

(3) For f,g£B, 

i d 3\^ (R ext (X±iO)f, g) 

is continuous. 

(4) If f is compactly supported, R ext {X ± i0)f satisfies the outgoing (for +) or 
incoming (for —) radiation condition. 

Proof. We prove the theorem for A + iO. We extend f £ B and u(z) = R ext {z)f 
to be outside Vt ext . Then it satisfies 

(H Q - z)u(z) = Ku{z) + j on Z d , 
where K = J2 n c n,P{n) is a finite sum of projections P(n) to the site n. Therefore 

(6.5) u{z) = R (z)Ku(z) + R (z)f. 

Let J be a compact set in Id, and take s > 1/2. We first show that there exists 
a constant C > such that 

(6.6) \\u(X + ie)\\^_ a < C||/||g, VA e J, Ve > 0. 

In fact, if this does not hold, there exists z M = A M + «e M , / M £ B, such that — 
Rext{Zfi)U satisfies 

(6.7) z^AeJ, II/mIIb^ ' = 1 as A 1 -> °°- 

One can then select a subsequence, which is denoted by {u^} again, such that 
converges weakly in TiT s . Since K is a finite dimensional operator, Ku^ converges 
in B. Therefore, in view of (6.5|, we see that converges in B* , hence in H~ s , to 
u such that ||u||^_ s = 1. It satisfies 

(H ext -X)u = 0, u = R (X + iO)Ku. 

Therefore u is an outgoing solution. By Theorem |6.2| u = 0, which is a contradic- 
tion. 
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We next prove that for s > 1/2 and f £ B, R ex t(X + ie)f converges strongly in 
V.~ s as e — > 0. To prove it, we consider a sequence = R ext {\ + ie^)f, e M — > 0. 
Then by the same arguments as above, one can show that any subsequence of {u^} 
contains a sub-subsequence {u^}, which converges in W.~ s to one and the same 
limit (independent of the choice of sub-subsequence) . This proves the convergence 
of R ex t(X + ie)f as e — > 0. Arguing similarly, one can also show that 

I d 3 \^R ext (\ + iO)feH- s 

is strongly continuous. The assertions of the theorem then follow from those for 
Rq(X + iO) and the formula 

R ex t(\ + t0) = Ro(X + i0)(l + KR ext (\ + i0)). □ 

6.2. Exterior and interior D-N maps. Let Hint — ~^disc + V be defined on 
Clint with Dirichlet boundary condition. The interior D-N map is defined by 



Ay (A)/ = d^Uint 



A £ a(Hint), 



where u int is the solution of the equation 
(6.9) (-A dtsc + V- X)u int = in 

The exterior D-N map is defined by 



int •> 



Hint 



/• 



(6.10) 



where 



A^(A)/ 



-d; 



A G Id, 



r(±) 



€ ;B* is the unique outgoing (for +) and incoming (for — ) solution of 



the equation 
(6.11) 



(H ext - \)u { ^] = in he 



v(±) 



on, 



/• 



The existence of u ext is shown by extending / to be zero on 7i d \dfl ext , and putting 

u e f t = / - R ext {\ ± *0)(-A dlsc - A)/. 

The uniqueness follows from Theorem 6.2. 

We represent u ext in terms of exterior and interior D-N maps. In the following, 
for a subset A in Z d , we use \A to mean either the operator of restriction 

(6.12) XA :^{Z d )3f^f\ 

A 

or the operator of extension 



(6.13) 



/, on A, 

0, on Z d \ A, 



which will not confuse our argument. We put 

C(i?o) = Oflint = 8Vl ex t = Clint n Cl ex ti 

1 d ^ 

( 6 - 14 ) •S'cfKo) = 4 E^(Ro)( 5 3 + (<SjT)XC(fl ). 

(Sju)(n) = u(n + ej ), {(Sj)*u){n) = u(n - ej ), 
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and also for n € C(i?o) 



(6.15) 



de sc(iio)( n ) = t *{ m e c ( R o) ; \ m - n \ = !}• 



For X e I d \ v{H int ), wc define the operator B^) R AX) £ B(£ 2 (C(Ra))) by 



(6.16) 4tk)W = A v( A ) - A S( A ) - A + ^gc^o) " ^c(Ko), 

where deg^^ is the operator of multiplication by deg^^^n). 

Lemma 6.4. Assume that X € Id\ o~(Hint), / € ^ 2 (C(i?o))- £ei ^ext an< ^ ^ e 
£/ie solutions of (6.11) and (6.9), respectively, and put 



X° 2 irit + X° «St + XC(R )f 



flint 



Then we have 

(6.17) ^)(n) = (i?(A±zO)xc(fl )4Uo)( A )/)( n )' n e Z "- 
In particular, 

(6.18) CtW - (£(A±iO)xc(fl„)5^ o) (A)/)W, n efiext, 

(6.19) f(n) = (R(X± l O) XC (R )B { c % o) (X)f)(n), n e C(ifc). 

Proof. Let r(n, m; A ± iO) be the resolvent kernel, i.e. 

r(n, to; A ± iO) = (fi(A ± i0)? m ) (n), 



where S m (n) = 5 mn . As in the proof of Theorem 5.8 by Green's formula, 
(6.20) 

53 ((A dlsc ^ ± ))(n)f(n, m; A ± iO) - u (±) (n)(A disc r)(n, to; A ± iO)) 

fie(ni„tua eit )nD(fl) 

= H ((^"'u^^WfCn.mjAiiO) - u (±) (n)(<9^«*f)(ri, m; A ± £0)) 



53 ((3""'!2 (±) )(n)f(fj, to; A ± iO) - (n)(d^ ext r)(n, m; A ± iO)) 

+ 5Z ((a^rf^^Hf^TOjAi^-^'M^rad^^TOjAiiO)) , 



n£D(R)\D(R-l) 



for sufficiently large integer R > 0. By the equations (6.9 1 and (6.11 ), the left-hand 
side of (6.201 is equal to 



(6.21) 



53 u^(n)((-A disc + V - X)r)(n, to; A ± iO) 

«£(ni„tun eal{ )nZ3(fl) 

53 u (±) (n)6„ m , 
ne(n int un ea! t)nr>(ii) 
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for any to G Z d . Note that, by our definitions of A^>(A) and h} ex ' t {X), 
9 n int ~(±) = d a int -. nt = A ~(A)/, 
= = -A^(A)/. 



The sum ^2 nGd Q int + Yl n edn m the right-hand side of (6.20) is then equal to 



n£C(R 



(6.22) 



53 ((A^(A)/)(n)r(n,m;A±iO) - f(n)(d^ int r)(n,m; A ± i0)^ 

" ) 

2 ((A^(A)/)(n)f(n, m; A ± iO) + /(n)(^?«*f)(n, m; A ± i0) 
:c(fl ) 

^ f(n,m;A±iO) Xc(Ko) (n) ((A?(A) - A<±>(A))/) (n) 

nGC(i?o) 

- 53 /(n) ((fl?*-' +fl?-*)F)(n,m; A ±t0). 



neC(flo) 

For n G C(Ro), 



((d^ int + d^)?) (n,m; A ± iO) 
- 53 m i A ± iO) — r(fc, to; A ± iO)) 



4 

o o 



(A^;5 C r)(n, m; A ± iO) — - 53 (r(n, to; A ± iO) — r(fc, to; A ± z0)) . 



4 

fceC(.R ),/c~n 



Therefore, the second term of the right-hand side of (6.221 is computed as follows: 



53 f(n>) (-A d i sc f)(n.,m;A±iO) 

neC(flo) 

+ \ E E (r(n,TO;A± l O)-f(fc,TO;A±zO) 

neC(flo) keC(Ra),k~n 

- - E T( n ) s n m + 51 ( ~ A + ^ de sc(fl ) w) /( n ) f ( n ' m ' A ± *°) 

neC(fl ) n6C(fl ) 



-J 53 r(k,m;\±iO) 53 /(n), 
where, in the 3rd line, we have used the fact that 



((- A dlsc - A)r) (n, to; A ± iO) = <5 nm , n, to G Z d , 
and exchanged the order of summation in the 4th line. Note that 



d 



E f(n) = J2((S j + (S j )*)xc( Ro )f)(k). 

n£C{R )-n~k j = l 
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Since we have for any m € D(R) 



E 



7(± 



(n)6 nm + ^ f{n)8 nm =u (± \m), 



o o 



neC(flo) 



(6.201 turns out to be 

S (± V)= E ^,m;A±iO)((A^(A)-Ai± ) (A))/)(r i ) 

neC(flo) 

+ E ( ~ A+ ^ de Sc(fi )( n ))/( n ) r(n,m; A±i0) 



1 d 
~4 E f ( n ' m ; A± * )E((^' + (^)*)^ c («o)/)H 

neC(R ) 3 = 1 



n€D(R)\D(R-l) 



for any m g D(R). In view of (6.16), we have thus arrived at 
u^(m) = (R(X±iO)xc( Ro )B^ Ro) (X)f)(m) 

+ E ( K {d rad u {±) ){n)r(n,m;X±iO)~u^ ± \n)(d rad r)(n,m;X±iQ)y 

n£D(R)\D(R-l) 

Taking the average of the sum with respect to R in the above equality, we have 
u^\m) = {R(X±iO) Xc{Ro) B i c ± { ) Ra) (X)f)(m) 



(6.23) 



+ ^ E {{drad-A ± {\,u n ))u^){n)?(n,m-X±iO) 



n£D(R) 



l - u ( - ± \n)({d rad -A ± (X,u} n ))?)(n,m- 1 X±iO) 7 



R 



ne-D(-R) 

up to a term of 0(R~ 1 ). By the radiation condition, we have 



1 

R 



i( d rad - A±(A, io n ))u (±) ) (n)r(n, m; A ± iO) 

neD(°R) 
/ 

^ E |((^ad-A ± (A,o; n ))n( ± ))(n) 



1/2 



< 



i? 

\ n£D(R) 



1/2 



— \r(n,m;X±- 



n£D(R) 



which tends to zero as R — > oo. The third term of the right-hand side of (6.23) is 
estimated similarly. This proves the lemma. □ 
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Lemma 6.5. Suppose A € /<* \ ir(ffmt)- Then for f,g£ £ 2 (C(Ro)), we have 
(6.24) (A?(A)/, d)e2 {C (R )) = (/. A v ( A )9)p ( c(i?. ))' 

( 6 - 25 ) ( A eS( A )/: ?) £ 2 (c(flo)) = (/, A e ^t ( A )ff)^ (c(fl[))) - 



(6.25 



(6.11 



Proof. The first equality (6.24) follows from Green's formula. We shall prove 
. Let u be the outgoing solution of (6.11), and v the incoming solution of 
with / replaced by g. For a sufficiently large integer R > 0, we have by 
Green's formula 

0= J2 (((-ffext - A)S)(n) • Wnj - u(n) ■ {{H ext - A)«)(n)) 

rie(_D(_R)nf2„ t )<= 



= ^ • «(n) + «(n) • (5? W £)(n)) 

nedD(R) 



As in the proof of Theorem |5.8[ we have 



^ [-(d^u)(n) ■ v(n) + u(n) ■ (d? {R) d)(n)) 

(-(dradu)(n) ■ v(n) + u{n) ■ (d rad v(n)j 



tedD(R) 



This implies 



neD(R)\D{R-l) 



= X! ((^ra<ju)(n) - A + (A,w„)u(n))w(n) 

T»£ZJ(°R)\D(fl-l) 



- ^ w(«)((9rod«)(n) - A-(X,u n )v(n)). 

neD(R)\D(R-l) 

Then, taking the average of the sum with respect to R, we have 
= ^ X! ((^rad - -4+(A,w„))u(n))w(n) 

n£-D(fl) 

- p X! "( n )(( 5 Tod - ^-(A,w„))u(n)), 

up to a term of 0(i? _1 ). By the radiation condition, we can see that the right-hand 
side tends to zero as R — > oo as in the estimate of (6.23). This proves (6.25). □ 
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7. Scattering amplitude and D-N maps 
7.1. Far-field pattern. We introduce the operator T^^A) by 
(7.1) f(±)(A) = G {± HX)xc^)B^ Ro) (X) : £ 2 (C(R a )) -+ L 2 {S d ~ x ). 

The main purpose of this subsection is to show that T^^A) is 1 to 1 (Lemma 

Although defined through ^/^(A), T^^A) does not depend on V. It is seen by 
the next lemma which follows from Lemma 16.41 and Theorem 14.71 



7.4) 



Lemma 7.1. Suppose A G Id \ o-(Hint)- Let be the solution of (6.11). Then 
we have 

+ o(|£r (d+1)/2 ) 

as I fcl — > oo. 



We need resolvent equations for R ext (X ± iO). Note that by (6.16) and Lemma 

OBaoW)* = Ay(A) - A^(A) - A + ^deg c(i? . o) - S c(flo) = fl^(A). 
Lemma 7.2. 

(7.2) R ext (X±iO) =R (\±iO)-R(\±iO)xc(R A%\\)xc(R )Ro(\±iO). 

(7.3) i? ea4 (A ± »0) = i?o(A ± iO) - i?o(A ± zO)xc(fl„)4^( A )xc(flo)^( A ± * )- 

Proof. Since u = R(X ± iO)x C ( i?o)j B^] ?o) (A)x c .( flo )_Ro(A ± iO)/ satisfies the 
equation 

(ff ea:t - X)v = in n ext , Volant = Ro(X ± iO)/, 
we have (7.2). Taking the adjoint, we obtain (7.3). □ 
We introduce the generalized Fourier transform in the exterior domain. We put 



F e xt (A) 



71 



: o(A) (l - Xc ( fl,)4Uo)( A )^(fl,)^(A ± iO)) , A e / d \ a p (ff mt ), 



and, in the same way as (4.18), we define 



Lemmas |4.6| and 7.2 imply that as |fc| — > co, 
(R ext (X±iO)f)(k) 

= e± (3-d)W4^| fcr (^l)/2 e±l /c., 00 (A^) a± ( A;Wfe)( g(±)( A )/)( ±Wfc ) 

+ o(|fc|-( d+1 )/ 2 ). 

This formula shows that G^t{X) does not depend on V . 

Lemma 7.3. For any <f> € L 2 (S d ~ l ), Qext{X)*4> satisfies the equation 



{H ext -X)G { -}{X)*^=Q m Q ext , (pUW**) 



an, 



= 0, 



and Gi x t(X)*(j) — Go(X)*(f> is outgoing. 
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Proof. By the definition, we have 



(A)> = (l - R{\ + ^0)xc(b. o )B { + { ) Ro) (X)xc( Ro) ) Go(X)*cb. 



By Lemma 



6.4 



v = R(X + iQ)xc\R„)B ( c { ) Ro) (\)xc(R )GoW* <t> satisfies the equation 

(H ext -X)v = in Vlext, v\an cxt = GoW*(f- 
The lemma then follows if we note that Go(X)*(j> satisfies 

(H ext - X)g o (X)*0 = Q in h ex t- □ 

Lemma 7.4. Suppose X E Id \ cr(Hi nt ). 

(1) f(±)(A) : £ 2 {C{R )) -> L 2 (S d ~ l ) is 1 to 1. 

(2) f(±)(A)* : L 2 ^- 1 ) -> £ 2 (C(i? )) is onto. 

Proof . Let us show (1). Suppose r^(A)/ = and let u^} be the solution 



that 



of (6.11). From Lemma 7.1 and the assumptio n, we have ~ 0. Then we see 



5.7 



By the unique continuation 



ju ext is compactly supported by Theorem 
property (see [IT], Theorem 2.3), we then obtain / = 0, which proves (1). This 
implies that the range of T^^A)* is dense. Since £ 2 (C(Rq)) is finite dimensional, 
(2) follows. □ 

7.2. Scattering amplitude. Recall that the scattering amplitude in the whole 
space is defined by (3.10). Passing to M\, we rewrite it as 

(7.4) A(\) = g( + \X)vg (\y. 

The scattering amplitude for the exterior domain is defined by 

(7.5) A ext (X) = ^ + HX)xc(r o) B { c \ Ro) (X)xc(r ) W)* ■ 
As in the case of Z d , we use its reparametrization on My. 

(7.6) Ae X t(X) = G i+ \X)xc(R a )B i c \ ) R JX)xc(R,)Go(Xr. 
Then we have as \k\ — > oo 

(Gi-h^rm)-(Go(xn)(k) 

(7.7) = -e^- d ^l^2i\k\- (rf-i)/a e <*-x 00 (A,o, i ) 0+(Aj u k ){A ext {X)4>){w k ) 

+ 0{\k\-^' 2 ). 



In fact, the left-hand side is e 



qual to -R(X + iO)xc(R a )B ( c !' ( ) Ro) (X)xc(R )Go(X)*(f). 
Using Theorem 4.7 we obtain (7.7). 

7.3. Single layer and double layer potentials. We have already introduced 
the operator i?(A ± iO)xc(R )BQ^ Ro ^{X), which is an analogue of the double layer 
potential. We also need a counter part for the single layer potential, which is an 
operator on £ 2 (C(Rq)) defined by 



C(R ) 



for / G £ 2 (C(R Q )). 
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The following lemma is a direct consequence of (|6.19|) and the fact that M^ R ^ (A) 



corresponds to Xc(Ro) R (* ± iO)xc(R )- 

Lemma 7.5. For A G Id\o~{Hi nt ), ^(^r )W^c(r )W * s identity operator 
on £ 2 (C(R )). 

7.4. S-matrix and interior D-N map. 
Theorem 7.6. For A G Id \ o-(Hi nt ), we have 
(7.8) A ext {X)-A{X)=f^\\)M^ ) {Xp-\x)*. 
As a consequence, S(X) and A^>(A) determine each other. 
Proof. Let us show (Q. For any (f G L^S"*" 1 ), let 

u = &-Xx)*4>-g£l{x)*4> 

= i?(A + i O)(xc(flo)^cUo)( A )^(i ? o) " 
In view of Lemma |7.3| u is the outgoing solution of the equation 

{Hext - X)u = in ttext, u\an Mt = G ( ^(X)*(j). 



(7.9) 



By (6.18), we can rewrite u as 



(7.10) 



u = R(X + *0)xc(«„)^ o) (A)xc(«„ ) a ( - , (A)> 



?(-)/ 



By (7.9), we have as \k\ — > oo 

u(fc) = C+|fcr( d - 1 )/ 2 e 2fe - a: ~ (A '" fc) a + (A,w fc ) 

x (^ (+) (A)(xc^„) J BS] ?o) (A)xc(fl„) - V)g o (A)*0)K) 
+ 0(|fc|-( d+1 )/ 2 ), 

where C+ = e ( - 3_!i - )7r2 / 4 v / 27r. On the other hand, by (7.10), we have as |fc| — ¥ oo 
u(k) = C+^-^-^^e^^^a+iX^k) 

x (G (+) (A)xc(flo) 4Uo) Wxctflo)^ (A)*?) K) 
+ 0(|fcr( d+1 )/ 2 ). 
These two expansions imply 

a (+) (A)(xc (% )<k)( A )^(«o) " V)G (Xy 
= Q {+) (A)xc(flo) Wxc ( r )Q { - ] (A)* . 

The left-hand side is equal to A ex t(X) — A(X). On the right-hand side, we insert 



»(+) 



after S^t^JA) to obt am 

= G (+ K^Xc( Ro) B^ Ro) (X)M^ Ro) (X) (B^ o) (A))*xc(H )^ ( - ) (A)* 
= f(+)(A)^ ) (A)f(-)(A)*. 
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We have thus proven (7.8). 



7.5 



Given A^(A), we can construct B^t^JX) by (6.16), hence M^ R ^(X) by Lemma 
Since T^^A) does not depend on V by Lemma 



A(X) b y Q . 

By (17. 81), we have 



7.1 



we can then construct 



f (+) (A)*(Ao: t (A)--4(A))f(-)(A) =f(+)(A)*f(+)(A)A4+] io) f(-)(A)*f(-)(A) 



Lemma 



7.4 



implies r( ± ^(A)*r( ± ^(A) is 1 to 1 on the finite dimensional space £ 2 (C(Rq)), 
hence bijective. Therefore, one can construct Mq? r JX) from A(X). □ 

8. Reconstruction from the D-N map 
In this section, we reconstruction V from the D-N map A^(A). 

8.1. Some properties of Schrodinger matrices. We identify — A<2i sc and A^(A) 

with matrices as follows. Let n( l \ ■ ■ ■ , are vertices in Qi nt and , • • • , 

are those in dVLi nt - We put 



■A/o = {: 



}, M = {,^ +1 \... , n H")}, 



and 



# {to G fi int ; to — n} = 2d, n£ fi int , 

44 ° 

^ # {to € flint ; m ~ n} = 1, n G 90 int . 

In view of the Laplacian on graphs, we construct a X (y+fi) matrix Ho = 

as follows (For the definition, see also [5]). 

1 



D = (dj. 



H = -(D-A), 

d eg ni „ t (^ w ) (i = j) 

1, if n« ~ for G fi in i or n^') G fi mt , 
0, if ?iW ^ n^, or ^),n^ G dto int . 



The potential is identified the diagonal matrix V = (vij) with 

.,.._/ ?(» (0 ) (i=j, *<v) 







(i ^ j or i > v + 1) 



Then = H + V corresponds to the symmetric matrix H = H + V. Moreover, 
identifying u with a vector (u(J\fo), u(A/i)) G C K+ ' i ) the equation 

(8.1) (-A disc + V)u = in fi int , 

is rewritten as 
(8.2) 

where by H(A/i; A/}) we mean a matrix of size #A/i x *A/}. The D-N map A^ 
is rewritten as 

(8.3) H(M;M)"(M) + H(A/i;A/oRA/o) = ?(^). 



HCAf ;M)«(M) + H(A/q;A/o)w(A/o) = 0, 



DISCRETE SCHRODINGER OPERATORS 



31 



Taking into account of the Dirichlet data 

(8-4) u\an mt = f, 

the above two equations are rewritten as 

<m> (« SKIS )-(&)• 

Assume that zero is not a Dirichlet eigenvalue of — A^ sc + V, which means that 



if u(Afi) = in (8.2), then u(Af ) = 0. Hence H(J\f ;AF ) is nonsingular. Then by 



using ( |8.2[ ), the D-N map corresponds to the fj, X fx matrix 

(8.6) A^/(M) := H(M;M)/(M) - H(M;M>)H(M); AA )- 1 H(AA ; M)/(M). 

To simplify the explanation, we translate O irlt so that 

(8.7) Cl int = {n e Z d ; 1 < n,- < M, j = 1, • • • , d} 
for a positive integer M. We put 

dnf = {nedn int ;n j = M+l}, dClJ = {n € dSU n t 5 «i = 0}, j = 1, • ■-,<*. 

Lemma 8.1. Given a partial Dirichlet data f on dCl int \ dClf and a partial Neu- 
mann data g on dd 1 , £/iere is a unique solution u on flint U 9f2j~ to the equation 

^ o 

(-A disc + V)u = in Clint, 
u = f on dCl int \dClf, 
tf}*>«u = g on df^. 

Proof. From the boundary values /(0, n 2 , ■ ■ ■ , H<j) and <?(0, n 2 , ■ ■ ■ , h<j), we can 
determine uniquely u(l, n 2 , • • • , n d ) for all 1 < rij < M for j = 2, • • • , d: 

u(l, n 2 , • • • , n<i) = -4#(0, ra 2 , • • • ,»<*) + / (°> n 2, • • • , 

From the equality ((— A^isc + V)u)(l, n.2, • • • , rid) = and the Dirichlet data f\gQ± 
for j — 2, - ■ ■ , d, we can compute u(2, n 2 , • • • , ftrf) as follows: 

1 1 d 1 - 

-2(2, n 2 , ■•• ,n d ) = -jzl "(i.^a,-- • ,nj + a,~- ,n d ) - -f(0,n 2 , ■■■ ,n d ) 

j=2 a=±l 

+ ^"(1, n 2 , • • • , rid) + V(l, n 2 , ■ ■ ■ , n d )u(l, n 2 , ■ ■ ■ ,n d ), 
for all 1 < nj < M, j = 2, • • • , d. We repeat this procedure to compute u(n) for all 

m = i,- •• ,m + l. □ 

For subsets A,Bc dCli nt , we denote the associated submatrix of A^ by A^(B; A). 

Corollary 8.2. Let u be the solution of (8.1), (8.4)- If f — on dCl int \ dClf , 
Ayf — on dCl^ , then u = in Clint- 

Corollary 8.3. The subm< 
dCl^ — > dfli is a bisection. 



V-> 

Corollary 8.3. The submatrix A^{dCt v ;dCl^) is nonsingular i.e. A^-(dCl^; dCl^) 
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Proof. Suppose / = on 90j n t \ dQ^ and A^/ = on dil 1 . By Corollary 



8.2 the solution u of 
implies that A^ (dQ^; dfl^) is nonsingular 



1), (8.4) vanishes identically. Hence / — on dQ^ . This 



□ 



Corollary 8.4. Given D-N map A^, partial Dirichlet data j% on d£li n t \ 90^ and 
partial Neumann data g on dfli , there exists a unique f on dfli n t such that f = fi 
on dVLi n t \ dili and A^-/| an - =g on <9f2^~. 

Proof. We seek / such that 



Av/| 



where h = f\ a n+- B Y Corollary 



i an; 
K = (A 



8.3 



v 



we take 

(dtv; dni))- 1 (g-A 9 (dni-,dn int \ dn+)f 2 



□ 



8.2. Reconstruction procedure from A^. We can now reconstruct V from A^. 
When d — 2, the procedure has been already given in [3], [4], [17]. For d > 3, we 
generalize this method as follows. 




Figure 7. The shape of Ci(0) in the case d = 3. 

We introduce the cone with vertex n € 51j„ t by 
(8.9) Ci(n) = jm e Clint ; ^ |m fe - n k \ < -(mi - ni)\. 

k^l 



If m satisfies the equation (8.8), we have 

(8.10) u{n) = c m u(m) 

meCi(n)\(n} 

for some constants c m . In particular, if u(m) = for all m £ Ci(?i) \ {n}, we see 
that u(n) = from ( |8.10[ ) (See also Figure [7]). 
Let II(p) be the rectangular domain defined by 

(8.11) n(p) = {(ni,--- ,n d ) e Q int ; m+n d =p, l<n»<M (2<»<d-l)}, 
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n i 



Figure 8. Situation of Lemma 153 



where M is from ( 8.7 1, and for r' — (r 2 ,--- , r^-i) £ [1, M] d 2 , we consider its 
section 

(8.12) I%r') = {{n^n' ,n d ) £ H(p) ; n' = /}. 

For d = 3, see Figure [8} 

Lemma 8.5. Assume M + 1 < p < 2M, and take a point (p - M - 1, r' , M + 1) e 
IT(p;r'). Lei u &e ifte solution of (8.8) with Dirichlet boundary data f such that 

( f(p-M-l,r',M+l) = l, 

\/ = on dSl int \(dntu{(p-M-iy,M+l)}), 
and Neumann data g = on dQ^ . Then we have 

!u(n) = if ni +w>d < Pi 
u{n) = if ni + n d = p. n' ^ r' , 
(8.14) u(p — M — 1 + i,r' , M + 1 — i) = (— 1)' /or p-M-l + i<M + l. 
If p = AI + I, taking the Dirichlet data f such that 

( T(oy,M) = i, 

\/ = on ffi jnj \(3O+U{(0,r',M)}), 
we /lave i/ie same assertion. 



(8.13) 



Proof. We put m = (p — M — 1, r', M + 1). First we show that m ^ Ci(n), if 
ri! + rirf < p. In fact, 

— (mi — ni) = n% — (p — M — 1) < p — — (p — M — 1) = M + 1 — n^, 
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and on the other hand, 

^ \m k - n k \ > \m d - n d \ = M + 1 - n d . 

k^l 

Then, in view of the condition for /, the Neumann data 9„u| an - 
we have u{n) = if n\ + nd < p. 

Assume that n\ + n c i = p and n' ^ r' . (See Figure [9j) Then 

— (mi — = M + 1 — nd- 

On the other hand, since n' ^ r' , we see that 

2J \m k - rifcl > |m d - rid | = M + 1 - n d . 

They imply m ^ Ci(n), hence u(n) = as above. 
Let us prove (8.14). Using the equation 



and (8.10), 



{(-A disc + V)u)(p — M — 1, r', M) = 0, 

and the fact that 

u = for m + nd < p, u(p — M — 1, r' , M + 1) = 1, 

we have u(p — M,r',M) = —1. Here we do not use the value of the potential 
V(p — M, r' , M). (See Figure 10 ) Repeating this procedure, we see u{p — M — 1 + 
i, r', M + 1 - i) = (-if inductively. □ 

Now we show the reconstruction procedure. 

1st step. We construct the boundary data / such that 

'/(M-l,r',M + l) = l, 

< / = on dn int \(dn+U{(M-l,r',M + l)}), 
A^f = on dD,^, 

y<j 



nr. 




Figure 9. Extension of the solution for the case (8.13). 
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by Corollary 8.4 Then the solution u of (8.1 ) and (8.4 1 satisfies the assumption of 



-1 (n = (M,r',M)), 

o 

(other n G Clint)- 



Lemma |8.5| By virtue of Lemma |8.8[ we have 

u(n) — 
Then, using the equality 



((-A disc + V)u){M,r',M) = 

and the boundary value f(M + l,r',M), we can compute the value V(M,r',M). 
Applying this procedure for all r', we recover V on all vertices (n\, r', rid) such that 
ni + n d = 2M. 

2nd step. Assume that we have recovered V on vertices such that n\ + n d > P for 
M + 1 < p < 2M. We construct the boundary data / such that 

f{p — M — l,r',M+ 1) = 1, 

/ = on dn int \(dntu{(p-M-i,r',M + i)}), 

A^- f = on dtti . 



By the same argument as in Step 1, the solution u of (8.1 1 and (8.4) satisfies (8.13), 



14). Since we have already recovered V on n\ + rid > p, we can compute u(ri) 
on m + rid > p using the equation (— Adi SC + V)u = and the boundary data /. 
Hence, using the equality 

((-A dtsc + V)u)(p -M-l + i, r',M + 1 - i)) = 0, 

and the fact that u(p — M — 1 + i, r', M + 1 — i) — (— l) 4 , we can compute V(p — 
M — 1 + i, 7'', M + 1 — i) for every i. Applying this procedure for all r', we recover 
V on all vertices {ri\, r', n d ) such that rii + n d = p with M + 1 < p < 2M. 
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3rd step. For p = M + 1, we construct the boundary data / such that 
' f(0,r',M) = l, 

< / = on dn int \(dQfu{(0,r',M)}), 
t A^ f = on dn^ . 



By the same argument as in Step 1, the solution u of (8.1) and (8.4) satisfies 

u(n) = { (n=(i,r',M + l-i)), 

\ (rii + < p or ni+n^—p, n'^r'). 

Then we can compute V(i, r', M + 1 — i) for every z as above. 
4th step. In the case n\ + n& < M + 1, we have only to rotate the whole domain. 
We have thus completed the proof of Theorem 
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